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Explicit  Formulae  for  the  Distribution  Function 
of  the  Sums  of  n unifomly  Distributed  Variables 

A„  fi5  Ostrovski # 


(1) 


1. 


Let 


«5  O <?  O $ 


be  n variables  subject  to  the  conditions 


(2)  ! y H,VO(v  - 

and  uniformly  distributed  In  the  corresponding  intervals.  %b& 
distribution  function  F(cr)  of  the  sms 


C 3 J ?X  t”  o o .+  § |2 

that  is  the  probability  for  the  inequality 


(*).  ^l  *..„+5n  i<r*  ^ 

s 

has  been  a as  far  as  I know,  given  explicitly  only  in  the  per- 

1) 

tieularly  simple  case  °<£  ®°^2  *.«3<x?j  by  Laplace  • Indeed  its 

expression  becomes /extremely  eates«  in  the  general  case  if 
one  does  not  use  a convenient  symbolism , 

2,  Me  us®  the  following  notations.  By  w©  denote 


Further  m will  use  the  "displacement  operator w &'  defined  by 

was  'prepared  under  contract  of  the  National  Bureau  of 
Standards  with  the  American  University,  Washington,  D.  C.  I 
should  like  to  gratefully  acknowledge  the  discussions  with  Mr. 

¥.  F.  Cahill,  Brs.J*  Ea  Curtiss,  ' C„  Elsenhart,  F;,  E.  Greenwood, 
A*  <J.  Hoffman,  E „ Litoes  and  J„  f odd . 

1 .Laplace ^collected  xmpers  vol.7,PP  257-260?cf  &1&®  JT.V .Uspensky, In- 
troduction to  Mathematical  Probabilities  1937, PP  277-278  end  the 
example  8-XO  on  page  200  jH  .Cramer  ,£-&thematleal  Methods  of  Statistics 
1940 ‘pp  2W-246*M.G.Eandftll*®ie  Advanced  Theory  of  Statistics, wol  1, 
third  edition , 19^7  ,pp . 140-145  sFeiler* Introduction  to  Probability 
Theory  and  its  Applications , vol  *X  ,Chap . XI , pp  .255  -236  9 problems  11, 


’v\.:£  \:o  , C ' ' -’v;^  i.u*  w4  ■•, 

\ ,.  ' - =••  ■ '.Vv  ; .. 


f 


m 


Sjj  f « s$  r Cm)  - r{n~ty . 


Then  our  expression*  for  Fjr)  is 


a.) 


(T)  F<«  - ,4p  - A ( > - X~  + < , S'*--  - -+  -■* ) 

3.  iTi  our  proof  of  (7)  we  eonslder^as  did  fcaplac©  (l.c)  5 
tho  corresponding  discrete  probability  problem 6 Lot  for  a certain 
integer  b isfeieh  will  toad  to  infinity^**  be  the  rational  fractions 


(8) 


j I f 4 m 


*h**e  the  integers  xv  are  "uniformly  distributed"  between  ~mv 
aad  ¥@  obtain  the  ezpr@gs.ions  Fa{sr)  of  the  probability  for 
the  S,tt@qm,l±tj 


(9) 


rt*- 


f rn  4 crm 


and  get -finally  Ffer)  &»  Xim  Fffi(o-)  if  a^-^oo  So  that 

(ID) 


si 


c< 


\>  *• 


4.  However  tafcan  the  above  notations  are  used*  the  standard 
method  of  complex  integration  gives,  as  ms  pointed  out  to  mo  by 


2*  The  corresponding  frequency  function  is 


I T * } F * (cr)  w ry; 


JX  (i-sia,d 4^) 


VI  *«  ^ 


1 * -%  ’•v® i v'  * K ’A  5 


this  follows  at  one®  from  the  relations 

(7W ) ^ ^ "t  < ~ n (<*  + 1 , ('  n - 3,  ~ ) 


a 

i?*»  ) f fa  + «r£*kr  - —V 

r t n * 5 


1 ^ ~ 1 /°s  , f ->* tv--) 


eh!  oh  v;t«  iaaedia  te  Xy  vs  rifled . 


-2- 


• ' 


. . 


B*.  Jo  If*  Curtiss,  a mrj  simple  proof  of  (7}  «, 

fj.  Another  approach  to  our  problem  which  goes  back  probably 
to  He  A*  Fisher,  19X5*  is  a,  geometric  one.  Denote  by  ? the 
parallelepiped  (2)  and  by  Rj»  the  common  part  of  P with  the  half 
plan©  (4)^ than  w®  have 


{11} 


M 

ITT 


where-  |Pl  « 2s o< . .c^^lsth®  volume  of  p and  j?r|  is  the  volume 
of  Pjy-  o Already  in  the  LapXacian  case  Kendall  (l.e)  has  used 
this  representation  of  p(<r)  and  derived, from  this  geometric  in- 
terpretatioa  another  proof  of  X&pX&ea8s  formula  in  using  a pretty 
involved  geometric  argument , It  ms  pointed  oat  to  sae  by  Dr, 

Au.  47  Hoffman  that  in  $sing  the  expression  (11)  the  formula  (T) 
can  he  easily  verified « and  Dr,  Hoffman's  proof  is  probably  the 
simplest  proof  of  the  formula  (7) * 

6*  However  it  seems  that  the  standard  statistical  method 
of  complex  integral  is  not  easily  applied  to  th©  following 
problem  of  a sore  special  character.  Consider  th®  a variables 
(1)  subject  to  the  conditions 

(12)  M&X  " 1.  £ 

then  we  mat'  to  determine  "the  probability"  F*(<r)  for  the  in- 
equality  (4) , However  th©  problem  1®  in  this  form  from  th©  point 
of  view  of  continuous  probabilities  not  even  completely  determined , 
It  is  'therefore  better  to  state  the  corresponding  problem  concerning 
discrete  probabilities , ¥e  consider  again  for  a positive  integer 
m a set  of  » variables  {8)  where  the  integer  numerators  rvra 
"uniformly  distributed"  between  -mv  and  mv«  Let  F%(ir)  be  the 
probability  of  the  Inequality  (9)  under  the  additional  conditio.? 


■ ■ ■ 

-•  -v» 


: . ■. 


■ 


■ - , •/  ■ 


(13) 


- 1 


Max 


then  F #{sr)  is  defined  as 


P*(cr)  - Lim  F%(oi 

fn  co  m 


W 

¥©  obtain  in  this  way 


(»v 


«o)® 


cis)  FVi.pfcnj  7X("*"‘'X«M 

^ °n  & v-.l 


r 

t 


To  This  expressless  ©an  be  interpreted  in  the  following 
m,j0  Let  ba  the  probability  of  the  inequality  (4)  under 

the  conditions  (12)  and  the  additional  condition  j 5^)«  • in 

assuming  that  and  are  equally  probable  ve  obtain 

fro®  the  formula  (7)  easily  ‘the  expression 

( ft)  rj(,  4 ^ _J1 

<16)  F ' (^)=  1 - S 

and  (15)  ean  b@  wit  tea  in  the  fora 


1 -tr  S ^ ;a  a*  \ / v ^ “ 4 

' — ^ A ( i + s ')(°<+v)t 

v = 1 


(17) 


P*(<r)  - A.°AA<r) 

r.=i 


(18)  % - -r_ ^£_  . 

-}-9. 

8&  It  sight  therefore  s©@®  that  (15)  can  foe  obtained 

by  straight  forward  application  of  the  elementary  probability 

theorems  If  each  co  is  interpreted  as  the  probability  for  the 

©quality  sign  5m  (2)  being  attained  just  for  -v  « ^ 0 However  the 

"weight  coefficients'5  co^  cannot  be  interpreted  in  this  way  since 

the  probability  that  ‘the  ©quality  sign  in  (5?)  is  assumed  at  least 

for  both  “v  =*  1 and  « 2 would  then  b©  4 <jx  lo  while  the  corres- 

^ ± z. 

ponding  terms  in  the  expression  of  F*{oi  obviously  tend  to  zero 

with  m ^ oo  0 


' 


V V;  ■■  ■ - ■ . ...  : , ; : • 


■ 


9*  On  t he  other  hand  the  geometric  interpretation  in  this 

mm  ls  sti11  possible , but  not  quit®  obvious*  ‘Th©  natural  im~ 

tarp^fetatiaa  would  fee  in  terms  of  the  (n~l)  dimsasXcm&X  surface 

_ , different 

■s^vds  os  f snd  Pg-..  However s/(n~l ) -dimensional  "faces*  of  these 

surfaces  would  have  to  fee  taken  with  the  weights , proportional  to 
their  areas. 


rf.g  a matter  01  fa,oty.t  was  the  problem  of  the  determination  of 
the  probabilities  F^fT)  which  arose  in  the  course  of  & discussion 
o.i.  a.  certain  modification  of  the  relaxation  method  and  was  the 


starting  po^fc  of  this  whole  investigation «,  For  our  ismaediate  purpose 
tfe®  p#C<r)  has  only  the  meaning  of  the  limiting  values  of  F*M{cr)!  , 

10,  It  my  further  be  mentioned  that  from  (xi)  it  follows 


iamaeaistely  that  F(cn  is  for  -«v#crJ-c<  a strictly  scmotonicaXlY 
i^SESSSiSfi.  function  without  stationary  intervals , same  then 


follows  for  the  expression  (16)  and  since  the  weights  (18}  are 
constants,,  we  see  that  F*(cr)  is  a strictly  monotonicaliy  increasing 
function  of  cr  for  - «*  & cr  4 ct. 


XI j Since  all  essential  points  of  our  original  proof  of  (J) 
occur  also  in  the  proof  of  (15) , 1 do.  not  give  explicitly  my 
proof  of  (?)  ? but  reproduce  in  section  12  with  Ms  kind  permission 
wi-v  proof  due  to  Br„  A*  J*  Hoffman  and  in  the  sections  13—22  give 
proof  ox  (15) » It  is?  believed  that  the  devices  used  in  this 


proof  can  be  employed  in  other  similar  .problems. 

12.  As  we  have  obviously  in  the  notation  of  no * 3 jp  j « £•%<%  . 


(?)  follows  at  once  if  we  only  show 

j_  ~fr  i 

n 5 * \ l i 


tt&) 


p i a 

!<r  * 


s*M-w 


v } ( Os  nr  cr 


>).. 


v'n 

! 

t 


' Ol  r-  Y , +• 


n / 


.a. 


V 


' 


f. 

■ ' t- 


, 


Fop  a » 1 the  right  hand  side  Is 


(20)  (a*+'c<i)+  - (c 

and  in  considering  separately  the  eases  CT^®<x,  Q&  rr  <<*, # 

~c<2.  c cr  < o and  cr^  we  see"  that  (20)  In  each,  ease  gives  the 
length  of  the  pant  of  C~,e<x^x>  with  x^rCF  * ¥©  can  therefore 
proceed  by  induction  and  assusa©  that  (19)  is  true  for  (n-1)  di- 
mensions • Consider  ’the  section-  PC'-4*)  of  with  the  (n«X)~dimansIomX. 
plan®  xn  **  u|it  is  obviously  gi/en  by  the  corresponding  coordinates 
xi  ? » which  satisfy  t 'ton  the  conditions 


(20) 


i & 
! ^ 


c* 


J 


v>-l 


^ Ot, 


n *1 


*,u+ + v,  cr-* 


therefore  ia  applying  (19)  to  p(u)  m have  for  its  {a-lj-dlxansiosal 

0” 

volume 


21) 


cr 


\ 

f h-i)t 


no 


S ‘ ’ v>)  ( ©C  - «<n  + o *—  U.  } 


n-1 


and  obtain  now 


(22) 


_J 

(nm  01 


Ao-s 


du. 


tout  the  integral  in  the  right  hand  expression  is  equal  toy  (7)  to 


* t [<*  ' "\n  O’-  0--«n  £j  ( I - ^ " }(f< + ^1 

13,  In  what  follows  we  will  have  to  use  the  "generalized 
toinoaial  coefficients M 

(23) 


n 

'-t  • 


y] . $ ( n / , ^ 0 

ni  v o y^e 


h - !j  a,  - - - a 


-6- 


-ws© 


. 

I. 


. 

, 


systematic  us©  of  this  notation  together  with  tie  sysaboli-sss 
introduced  in  130  „ 2 appears  to  be  very  convenient  In  many  coaabina-* 
torial  discussions.  Its  use  would  also  simplify  considerably 
Laplace fs  argument  (l.c.)  for  the  case <=<,  esc*^  |>re  fe&ve 

fcr-  any  integer  !c  ^2  and  any  integer  Y 


!s,)  l *H  dj'-  [d]  ^ 


; 


as  • follows  at  once  from  {23)  for  and  f-ko. 

In  applying  (24)  repeatedly  we  have  for  two  integers  kla  k2* 
and  any  integer  K i 2 


(»)  £ UJ-  t\"]-  [ J].  (1-5 

Vsh«  ^ h J 


<.  **.-  Ki  'tli  r 


] 


14 » Suppose  now  that  we  have  for  at-»eo  and  a variable  in- 


feeger  to 

(25)  - 

| -»(3  '•  ^>00)  Sj 

then 

(at) 

la  j/lr-  -* (3^!  (n  -1,2,...)  . 

Indeed  if 

QkQs  then  | ^ j is  either  aero  of 

{28} 

.1  1 » 
n!  s ks 

And  if  (5  > 

Oj,  the  left  side  in  (28)  tends  to  77  6 Idri  generally 

if  yo  assume  In  addition  to  (26)  that/a  variable  Integer  X r ^/m 
then  370  have  from  (27) 


-7- 


> 


(29) 


--“■4 


n A-  ...  n 
6 - X / , - 0 pj 


't 


+ 


m 


13 


t K 


5*  ■ We  have  with  the  notation  (2 


(?o$  (i-x)-ft  - T 

p2-'®0 


f>4B~X  1 

U»“l  J 


A? 


\ h ^ **  u & o ^ ° J 


Indeed  the  right  side  expression  is  by  (23) 

co  -1 


n-2. 

rl— IX 


"»Ji  © 


4 X- 


\si”  S 


(31) 


® (X-x) 

3.6  u If  we  have  a development 

t?o 

?(x)  «*  21  r(h)xk 


K*-«® 


we  have  fox*  say  integer 


OO 


xyPc«)’  21  w<K)x 


K +•  'I 


oo 


? NU-f)*^  2L  S w (*)  * 

-fionam  i i 


*4 


,3.  = - «> 


fv-  - OO 


and  more  generally  for  any  polynomial  in  z 
(32)  #(x)P(x)  -Z.  J#(Sk  )H(* )xK  . 

K* 

17*  We  consider  now  a integers  (n  \2)  *i8p* ,. . , ,?%  and 

denote  by  Ok  the  number  of  representations  of  t he  Integer  k as 
the  sine 


( 33)  (Ofe)  k ® xx  4* . * o4*xn  (jxjj  | • « » .*  1 


x. 


n 


1 & 


of  ii  integers  x-v  satisfying  the  indicated  conditions,  Then  w©  can 


write  :0t  « mn  4 „ . .+m^2 

cj© 


2l  t *s- 

•K*-«  •' 


• inrV 


- m 


* s 


{#  4 X 4-«  - + 

ar»v4  .1 


X “ + X v „ _ _+■  X / - -. 


■"V  * 


A'") 


x 


f'\  ( i+  + x: 

Vsi 


m*>  - x~M  A (-7-iT-  ) 

' V=1  * * ? 


•8- 


' J. 


' 


(3*0 


O© 

c«x  “ 

;%  - oc 


y-M,  Nn  JL  , 

x (i-x)  7T  ( >- X v /; 


the  right  eld©  expression  is  by  {30}  ma  (32) 

V - J.  OO 


n 

7T 


XiBvH 


00 


F) 


= E 5K  7T 

then  ws  ofci^n  "*  ‘v<*i 


a mv*l  -V  \ <\  -*  Vi-ll  k 

S * ? I X ® 

°K  ' L «-*  J " 5 


(35)  |C  = 

* V**  H 


fnv  + i>  [ ^ -ir  W.-S-  ^ i* 


s.  0 - i 


1 «IL  , 3.  m v*  * \\  lv\*  V\--  H -»i  j 

i = A d r>-V  j 


l8o  From  (35)  we  have  again  for  two  integers  ^ y ^4 4 


L. 

and  since  the  right  hand  sm  1©  by  (25) 


ic  =I  7t(*-s^r,)[^;1ni=  A 0-s*""1") 

k*K«  ' -Jai  ■'  * J ">>*1 


, K*  r . 

I m v>  e \ : HtMn-1 


/ 


K a ft  4 


h - 1 


r 

! 


?i>Ki+  n 


r 


+ {(^4  “h 

n 


- ( 


K ~ «■  V\ , ^ i.  \ C + K 1 1 1 

/ 1 n 


grs.i‘:\ 


Z C 

ii\“  &.S. 


s rvl)7\( 


a TA  v -t  1 


M 


H-f  n 


19.  ¥©  denote  now  by  0 the  number  of  representations  of  k 


as  the  sue  of  n integers  *?  satisfying  the  conditions  )xv|<  m 


f 


V «*  15  toe*} 


Obviously  we  obtain  S 3&  from  Cjr  ^L*VA  replacing  each 


<V  ■ 


' ■ 


\ ■ : . 


<W 


X.f  m denote  now  by  % the  number  of  representations  of  k as  the 
«?  a Integers  x^  satisfying  th©  condition 

(;>8)  Max  ^ ,1 

V m v 


^ obviously  Die  «•  C2c  - and  from  (35)  and  iTT) 


20.  In  order  to  obtain  now  the  probability' ?*ffi(9-)  of  no.  6 
w obtain  th®  numbar  of  all  "favorable  oases”  in  replacing  la  th© 
right-  head  expression  of  (39)  by  an  arbitrary  negative  cumber 
# and  kg  by  the  greater  Integer  k « for  mj  contained,  in  srjae  Then 
we  obtain 


<*>)  Nm(<rt*.  J f\  (i-S**’’4)-  S"t\  (i-s'""'1  ■1)?[ 

l r‘  J 


®»  of  "all  possible  oases”  Is  obviously 

*>  >n 

im  * 7T  ( *.  •'»„  i-i)  - 7T  ( *mv-i'5  % 

^ 3 § Vs  1 

by  dividing  we  obtain 


$v  „A 

C4a,t  ^(g-;*=  ~§; 


- ~7T 


w 7\  C a w.  v+i)-  f \ (im~  i) 


S 43L  am^lw  n-S-y 

|A('-Sr  • ^AK  4 n 


= 1 


V=i 


obtained  now  the  solution  of  the  discrete  probability  problem 
fo:rnuXated  in  no.  6. 

21 » She  expression  !sf  (or)  can  be  brought  into  a form 
'it&jth  ”ssefuX  as  well  for  a geometric  interpretation  aa  for  the 
passage  to  the  limit . Wo  can  obviously  write  for  the  expression 
:»wm^  ‘the  braces  ia  (40) 
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fo  obtain  now  the  limit  of  the  expression  (42)  m divide  both 
%(tr)  and  & ffi  by  For  the  denominator  w©  obtain 
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ih:b?»  transformation  la  obtained  in  special i sing /?•  h.?  appropriate 
way  the  Ida  ntity 
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Is  developing  the  term  corresponding  to  a value  ©f  ^ w;  obtain 
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where  a is  an  Integer ^ -a  <&  <n.TSae  general  term  of  tim  right  hand 


side  converges  fcy  the  formula®  (29)  and  (46)  to 
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and  the  limit  of  (49)  1®  therefore 
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She  formula  (15)  follows-  Immediately, 
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portion  of  the  NBS  work  is  performed  for  other  government  agencies, 
particularly  the  Department  of  Defense  and  the  Atomic  Energy  Com- 
mission. The  functions  of  the  National  Bureau  of  Standards  are  set 
forth  in  the  Act  of  Congress,  March  3,  1901,  as  amended  by  Congress 
in  Public  Law  619,  1950.  The  scope  of  activities  is  suggested  in 
the  listing  of  divisions  and  sections  on  the  inside  of  the  front 
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from  the  Government  Printing  Office : the  Journal  of  Research,  which 
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Technical  News  Bulletin,  which  presents  summary  and  preliminary  re- 
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($1.00).  Information  on  calibration  services  and  fees  can  be  found 
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(25  cents).  Both  are  available  from  the  Government  Printing  Of fice. 
Inquiries  regarding  the  Bureau’s  reports  and  publications  should  be 
addressed  to  the  Office  of  Scientific  Publications,  National  Bureau 
of  Standards,  Washington  25,  D.  C. 
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